THE GENERALIZED STAR PRODUCT AND THE FACTORIZATION OF 
SCATTERING MATRICES ON GRAPHS 
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Abstract. In this article we continue our analysis of Schrodinger operators on arbitrary graphs 
given as certain Laplace operators. In the present paper we give the proof of the composition 
rule for the scattering matrices. This composition rule gives the scattering matrix of a graph as 
a generalized star product of the scattering matrices corresponding to its subgraphs. We perform 
a detailed analysis of the generalized star product for arbitrary unitary matrices. The relation to 
the theory of transfer matrices is also discussed. 



1. Introduction 



Potential scattering for one particle Sciirodinger operators on the line possesses a remark- 
able property concerning its (on-shell) scattering matrix given as a 2 x 2 matrix- valued function 
of the energy. Let the potential V be given as the sum of two potentials Vi and V2 with dis- 
joint support. Then the scattering matrix for V at a given energy is obtained from the two 
scattering matrices for V\ and V2 at the same energy by a certain non-linear, noncommuta- 
tive but associative composition rule. This fact has been discovered independently by several 
authors (see e.g. [H], |7[ ||, g |65|]) and is an easy consequence of the multiplica- 

tive property of the transfer matrix of the Schrodinger equation (see e.g. pi]]). It has also 
been well known in the theory of mesoscopic systems and multichannel conductors (see e.g. 
[71, 22, 23, 24, 25, 55, 0, EO, 28|). In higher space dimensions a similar rule is not known. 



However, for large separation between the supports of the potentials the scattering matrix at a 
given energy is asymptotically related to the scattering matrices for Vi and V2 at the same energy 



[39, 40] 



To the best of our knowledge the composition rule for 2 x 2 scattering matrices was first 
observed in network theory by Redheffer [|^, 58 1, who called it the star product. In our 
preceding article [ p2| ] we extended this result to quasi-one dimensional quantum systems - 
Schrodinger operators on graphs. Such systems are nowadays a subject of intensive study (see 
e.g. [ |3T| , 0, 13, 14|). Some other related works are quoted in [42]. In [|^ 



48 



17, g]. Scattering 



, |67p differ- 
ential operators with Neuman boundary conditions on "fat graphs" were considered, i.e. on thin 
domains in M'^ which asymptotically shrink to a graph. 

There is also a large amount of literature on linear difference operators on graphs. The moti- 
vation for the study of such operators comes from the graph theory, where the spectrum of these 
operators are known to be related to topological properties of the graph [ 
theory for such operators was developed in [ p6[ , ^. 

In [^] we considered the (continuous) Laplace operator on graphs with an arbitrary number 
n of open ends (i.e. channels) and with arbitrary boundary conditions at the edges resulting 
in a self-adjoint operator. We formulated and proved necessary and sufficient conditions for 
such operators to be self-adjoint. We provided an explicit expression for the resulting unitary 
nXn scattering matrix in terms of the boundary conditions, the lengths of the internal lines and 
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the given energy. Furthermore, we generaUzed Redheffer's star product to what we called the 
generahzed star product. This is a non-linear, noncommutative but associative composition rule 
for unitary matrices not necessarily of equal rank and resulting in a unitary matrix. 

Under special circumstances there is an alternative way to describe the generalized star prod- 
uct. Fix p>l and consider the group U{p, p) with its natural multiplication. As a set this group 
is isomorphic to some subgroup of \}{2p). This non-linear set isomorphism is well known in 
the case p=\ (see e.g. [26]) and can be easily generalized to the case of arbitrary p>\. Under 
this isomorphism the multiplication in \}{p,p) induces new nonlinear multiplication *p in this 
subgroup of U(2p), which is our generalized star product. The operation *p can be extended by 
continuity to the whole U{2p). The set \}{2p) with *p as multiplication is no longer a group, but 
only a semigroup. 

Employing this generalized star product in [ p2| ] we provided a formal proof based on the 
quantum mechanical superposition principle to show how the scattering matrix at the same en- 
ergy for the whole graph can be obtained from the scattering matrices of two subgraphs obtained 
by cutting the graph in any way in two. Again for the special case of 2-channel scattering ma- 
trices, like potential scattering on the line, this formal argument is well known (see e.g. [20]). 
In this article we will provide a rigorous proof of this composition rule. It is interesting to note 
that in this general case the composition rule cannot be reduced to the multiplicative property of 
the transfer matrix of the Schrodinger equation on the graph. 

Such composition rules are important in the study of the electric conduction in multi-terminal 
mesoscopic systems. By the Landauer-Biittrker theory the electric conduction in mesoscopic 
systems is directly related to the transmission probability and thus to the scattering matrix [ ^1] , 
p^ , pl| , 12]. A good introduction into the theory of electronic transport in such systems is given 
in the book [ pO| ] by S. Datta. The formal arguments leading to the composition rule for the 
scattering matrices are presented on p. 125 - 126 of this book. 

The composition rules are also very useful in the study of statistical properties of large random 
or periodic systems. Examples of such systems can be found e.g. in [|, 0, |27|]. In [|l|, || 



44 , 46] we proved that in arbitrary dimensions the scattering phase (or more generally of the 
spectral shift function) per interaction volume equals (up to a factor %) the difference of the 
integrated densities of states for the free and interaction theories respectively. In the strictly 
one-dimensional situation (Schrodinger operators on the line) the Lyapunov exponent is known 
to be related to the logarithmic density of transmission probability [52, 53, 54, ^]. Due to 
the Ishii-Pastur-Kotani theorem (see e.g. [p^) the vanishing of the transmission amplitude for 
almost all values of energy implies localization (i.e. the spectrum must be purely point), see also 
the related works [0, || |2|, |4[ || ||, mi]. 

Certain Laplace operators on (infinite) periodic graphs were previously considered in [}7|, |27|]. 
There are also some attempts to consider differential operators on regular graphs with random 
boundary conditions or on random graphs with deterministic boundary conditions (see e.g. [^). 
Some other examples can be found also in [38, Chapter 3]. A difference Laplace operator on 
the edges of aperiodic tilings was considered in [37]. Such systems provide a main field of 
application for our composition rule which will be discussed in a forthcoming publication. 

The article is organized as follows. In Section ^ we recall the general quantum scattering 
theory on graphs as given in In Section I we recall the definition of the generalized star 
product and study some its properties. In particular we show that this product applies to arbitrary 
unitary matrices. In Section ^ we give a rigorous proof of the composition rule for scattering 
matrices on arbitrary finite graphs. Section ^ is devoted to the special case of graphs having an 
even number 2p of external lines. If the new graph is obtained by gluing of exactly p lines then 
it has again 2p external lines. We consider the question whether in this case the composition 
rule for the scattering matrices can be reduced to the multiplication rule of the corresponding 
transfer matrices. In general for p> \ the answer is negative. We formulate a necessary and 
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sufficient condition, which guarantees that the composition rule for the scattering matrices is 
equivalent to the standard multiplication of transfer matrices. 



After completing the work we received the preprint [34] by M. Harmer where among other 



questions the composition rule for the scattering matrices is also considered. The results there 



partially recover our Theorem 4.1 below 



We are indebted to P. Kuchment for sending us the preliminary version of the preprint [48] 



and also for pointing out the works of R. Carlson 1 13, 14, 15]. 

2. The Laplacian on a Graph and its Scattering Matrix 
In this section we will recall the definition of Schrodinger operators on an arbitrary but finite 



graph and the construction of their scattering matrices [142[]. 

We consider an arbitrary graph F with a finite number « > 1 of external and a finite number 
m>0 of internal lines (edges). More precisely this means that outside of a finite domain the graph 
is isomorphic to the union of n positive half-lines. Any internal line ends at two, not necessary 
different vertices and has a finite length. We assume that any vertex of T has non-zero degree, 
i.e. for any vertex there is at least one edge (internal or external) with which it is incident. 

Let the set £ label the external and the set J the internal lines of the graph. We assume that 
the sets £ and J are ordered in an arbitrary but fixed way. To each e g £ we associate the infinite 
interval [0, oo) and to each /g J the finite directed interval [0, a,], where a,>0 is the length of this 
line. With this association the graph becomes directed, such that the initial vertex of an edge of 
length a corresponds to x = and the terminal vertex corresponds tox = a. The external lines are 
assumed to be directed in the positive direction of half-lines. 

We define the Hilbert space "K = L^iT) as 

ee £ ie 3 

where IKg = L?{0, °°) and [K, = L^(0, a,). Elements of !H are written as column vectors 

(2.1) V= ({\|/Jes£, {¥/}fe3)^ = (¥£'¥3)^' ^e^^e, \|/,- G IK/. 

Similarly we define the Sobolev space IV^'^(r) as 

W^'^iV) = W^-^iO, oo) ® W^-^iO, at), 

ee £ ie 3 

where W^'^{0, oo) and W^-^{0, a,) are the usual Sobolev spaces of square integrable functions 
whose distributional second derivatives are also square integrable (see e.g. [|9|]). Let [ ] : W^'^(r) - 
(^2(n+2m) surjcctivc linear map which associates to each \|/ the element [\|/] given as 

again viewed as a column vector with the same ordering as in \|/, i.e. with the ordering given by 
the ordering of £ and 3. 



In [42] we showed that for any two (« + 2m) x{n + 2m) complex matrices A and B with AB* 
being Hermitian and the {n + 2m) x 2{n + 2m) matrix (A, B) having maximal rank equal to « + 2m, 
one can define the self-adjoint Laplace operator A(A, B,q) in "K corresponding to the boundary 
condition 

(2.2) A\|/ + B\|/' = 0. 

Here £=(£?!,.. . ,<3„,)^g ]R™,m = #(J). Furthermore, any self-adjoint extension of the Laplacian 
on the given graph is given by A(A,B,a) with some matrices A and B satisfying the properties 
stated above. If J = we simply write A(A, B) instead of A(A, B, ■). 

Before we turn to the scattering theory for A(A, B, a) we recall some well-known facts from 
scattering theory in two Hilbert spaces and S)2 (see e.g. []60|]). Let Hi and H2 be self-adjoint 
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operators in the Hilbert spaces and ^2 respectively. Let J be a bounded operator from ^1 
into Sj2- The two-space wave operators are defined as the strong limits 

Q.HH2, Hu3) = s-hm e'"'-'de-'"''PAHi ), 

where PudH) denotes the projection onto the absolute continuous subspace of H. We consider 
the sets 91+ of elements g & S)2 for which limj^ipoo \\d* e~'^^' P^c{H2)g\\ = 0. The wave operators 
Or are called 0-complete if i^z = Ran(n*)®OT+. For details we refer e.g. to Chapter 3 of the 



book [72] 



Now we consider a graph consisting of the external lines of the original graph F. On 
the graph Fg we consider the operator -A{A^ = 0, = I) corresponding to Neumann boundary 
conditions. Let "3: Ji ^ JCg be given as = \|/£ according to the notation (^). In partic- 
ular is identity if m = 0. Since we actually deal with finite dimensional perturbations by the 
Kato-Rosenblum theorem (see e.g. Theorem 6.2.3 and Corollary 6.2.4]) the wave operators 
n-(-A(A, B, a), -A(A£ = 0, = I); J) exist and are ^-complete. Thus the scattering operator 

(2.3) S(-A(A, B, a), -A{A^ = 0, = I); d) = {Dr)*Q.+: IKg ^ Ji^ 

is unitary and its layers SA,B.a{'^)- C" —>■ C" (in the direct integral representation with respect to 
-A(A£ = 0, Bg = I)) are also unitary for almost all energies X g R+. 

The resulting scattering matrix is related to the scattering wave function for the operator 
-A(A, B, a) at energy X>0 as follows. The function \|/*^(-, X) indexed by k& £ and with components 



(2.4) w){x,X)=< 



Sjk{X)e''^'' for j^ej^k, 

g-,VIv + Skk{X)e'^' for j G £, j = it , 

^ ajk{X)e'^' + ^jk{X)e-'^'' for j g 3, 



solves the Schrodinger equation with the operator -A(A, B, a) for energy X>0. 

Recall that in potential scattering for one particle Schrodinger operators the wave operators 
give precise meaning to the scattering wave functions, i.e. solutions of the Schrodinger equation 
at positive energy X>0. Similarly the wave operators D.-{-A{A, B, a), -A(A£ =0, Be =1); d) deter- 
mine the "external part" of the scattering wave function, i.e. \|/^(x, X) for j g £. The completness 
of usual wave operators means that any solution of the Schrodinger equation at energy X>0 can 
be uniquely represented as a superposition of the scattering wave functions. Similarly, in the 
present context the ^J-completness of the wave operators D.-{-A{A,B,a),-A{A^ = 0, Bg = I);S) 
means that the external part of any solution for the Schrodinger equation with the operator 
-A(A, B, a) at energy X>0 can uniquely be represented as a linear combination of the external 



parts of the scattering wave functions (2.4) 



The scattering matrix SA,B,aO^) as well as the mxn matrix amplitudes aA,B,fl(^) and ^A,B,d^) 
are. determined as solutions to the equation 

fS{X)\ ft 

(2.5) ZamA'^) a(X) =-(A - iVXB) 

where 

(2.6) ZA,B,aCX) = AXaiX) + iVXBYa{X) 
is the {n + 2m) x{n + 2m) matrix with 
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and 



I 

Here g-'^^ stands for the mxm diagonal matrix with elements 

If det ZA,B,aO^) ^ the scattering matrix S{X) = Sa,bm(^) well as the mxn matrices a(^) = 
ctA,B,a(^) and = (3/i,B,a(^) can be uniquely determined by solving the equation ( ^^ ) in the 
form 

fS{X)\ /I' 

(2.7) a{X) = -ZA,B,a0^r^ (A - iVXB) 

We denote by T.A,B.a = > | det Za^b,JX) = 0} the set of exceptional points for which Za^^JX) 
is not invertible. 

Let (|) be an arbitrary measurable function on the graph T. The subset supp (j) of all edges of 
the graph T will be called the support of the function (|) if (|) ?i a.e. on supp (j). 
In p^ ] we proved the following 

Theorem 2.1. For any boundary condition {A,B) and arbitrary a g the set '^a.b.u equals 
the set OA,B,a of all positive eigenvalues of — A(A, B, a). This set is discrete and has no finite 
accumulation points in R+. The eigenfunctions of—A{A,B,a) have support on the internal lines 
of the graph. For all X g '^A,B,a the equation CLS) is still solvable and determines SA,B,a(X) 
uniquely. 

Given an arbitrary nxn unitary matrix U and an energy ?io>0 we can find boundary conditions 
A,B defining a self-adjoint operator -\{A,B) on a single-vertex graph (i.e. with m = 0) with n 
external lines such that the corresponding scattering matrix is given as SamC^)) = ^- The proof 
of this fact can be found in [45]. For other inverse problems on graphs we refer to [32, 15]. 

Recall that by definition the operator A{A,B,a) is real if it commutes with complex conjuga- 
tion, i.e. for any \|/g D{A{A, B, a)) the function \|7 also belongs to D{A{A, B, a)) and A(A, B, a)\|/= 
A(A, B)\|7. Equivalendy, this means that any \|/ g D(A(A, B, a)) (i.e. \\f g W^-^iF) satisfying A\|/ + 
= 0) also satisfies the equation A\|/+B\|/' =0. Thus, A(A, B, a) is real iff Ker(A, B) = Ker(A, B). 
The last condition is satisfied iff there is an invertible matrix Ci such that A = CiA, B = C\B or 
alternatively there is an invertible matrix C2 such that both C2A and C2B are real. We recall that 
A(A, B, a) = A(CA, CB, a) for any invertible C (see [g2|]). 

In [^, Corollary 3.2] we have proved the following 

Theorem 2.2. For arbitrary a g W", X>0, and all boundary conditions A, B defining the self- 
adjoint operator A(A, B, q) 

(2.8) SjsjXf = SA,B,aa). 

In particular, if the operator A(A, B, a) is real then SA,B,a(X) = SA,B,a(X)^ for all X>0. 

Here we give an alternative proof. 

Proof. From the selfadjointness of the operator A(A, B, a) it follows that the matrix A + iVXB is 
invertible for all X > (see [42]). The relation ( ^ ) implies that 

Sa.bA'^) \ I I 

o^A,B,a(X) = -(A + iy/XB)~^ (A - iy/XB) ^A,B,a(X) 
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Similarly, for the operator A(A, B, d) we have 



a- 



lB,A) = -(A + iVXB)-^ (A - iVXB) pA,fi,«(^) 



a- 



Taking the complex conjugate and multiplying by (A + iVXB) '(A - iV^B) from the left we 
obtain 



Pa = -(A + iVXB)-^ (A - iVXB) 



-i\f\iif^_ 



We multiply this relation by ^^^(X)^ from the right and make use of the unitarity of the scatter- 
ing matrix in the form S-^-g{X)S-^-g(kY = I thus obtaining 



^A,B,a 



/ 



= -(A + iVlB)'^ (A - iVlB) 



V^^^^)Sa:bM 



Equivalently, this relation can be written in a form analogous to (2.5), 



(2.9) 



ZA,B,a(^) 



/ ^a:b,c0^ 

\^A,B,a(^)^A,B,g(^)^ 



= -{A-i 



ft 

iVlB) 



In [ p^ ] we proved that the equation ( |2.5| ) has a solution for all X > 0. If X > is not an eigenvalue 
of the operator A(A, B, a) then it has a unique solution. If X > is an eigenvalue of A(A, B, a) 
this solution is non-unique, but still determines the scattering matrix uniquely. Therefore from 



comparison of ( |2.5[ ) and ( |2.9| ) the relation ( [2.8D follows. If A(A,B, a) is real by the remark 
preceding the theorem we can choose the matrices A and B to be real. From this and from (2.8) 
the second claim of the theorem follows. □ 



We note that the comparison of ( |2.5[ ) with ( [2.9| ) also gives the relations 

CA,B,a(^) = PA,B,a(^) SA,B,aQ^V < 



pA,B,fl(?l) = aj-gjl) SA,BA'>^f ■ 

3. The Generalized Star Product 

Let U^^^ and ?7^^^ be arbitrary ni xni and ?i2 unitary matrices respectively. Let p be some 
integer satisfying l< p<{ni+ n2)/2, p < nj, j =1,2, and V be an arbitrary px p unitary matrix. 
We write U^^^ and U^^^ in a 2 x 2-block form 



(3.1) 



U 



(i)_ 





where and ?7{j' are pxp matrices, u[\^ is an {ni-p)x{ni-p) matrix, is an {n2-p)x(n2-p) 



matrix etc. The unitarity condition for ?7*^^ then reads 



(ir 



11 u\ 

(ly 
12 



(1) 

1 ^"^21 '^21 



'22 



(1) 

'22 



Ujj +'^21 "^22 

u^2 ^^22 ^21 



I, 
I, 
0, 
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and similarly for U^^\ 

Definition 3.1. The unitary matrix U'-^^ is called V-compatible with the unitary matrix U^'^^ if 
the py- p matrix y?722'^*^i? does not have 1 as an eigenvalue. For the case V = 1 the matrix 
U^^^ is simply called compatible with U^^^ (for given p>l). 

Note that the compatibility of the matrices is not a symmetric relation, i.e. if U'-^^ is V- 
compatible with U^^^ then U^^^ need not be V-compatible with 

Obviously, if U'-^^ is V-compatible with U^^^ then also the matrix V*u[^iVU22 does not have 
1 as an eigenvalue. Indeed, let us assume the converse, i.e. let there be non-zero c e such 
that V*ull^VU^2^c = c and thus 

VU^^V*uiyU^^c = VU^^c. 

Since VU22C'^0 the matrix V?722'^*^if has 1 as an eigenvalue, which is a contradiction. From 
this it also follows that if U^^^ is not V-compatible with U^^^ then the matrix V*u[^VU22 has 1 
as an eigenvalue. 

From the unitarity conditions it follows that 

o<u[fu[\' = i-u^fui\'<i, 

0<<^<^I-«<I, 

and thus ||t/i7||<l, ||t/22^||<l. Similar inequalities hold for and Therefore ||V[/2T^*^i?ll^ 
1. Strict inequality holds whenever ||?/22^|| < 1 or < 1 and then U^^^ is V-compatible with 

U^'^^ and U^^^ is V-compatible with ?7*'' for all unitary px p matrices V. 

In general if U^^^ is V-compatible with U^^^ then it is easy to see that the following p x p 
matrices exist: 

^1 = (I- v[/^2^v*[/*fr'v= v(i-[/^2'^*[/{iVr\ 
^^'^^ K2= (i-v*[/{fvc/^2'riv*=v*(i-[/{fv[/^2'^*ri. 

An easy calculation establishes the following relations 

Ki = V + VUl2^V*ulfKi = V + VU^22K2Ul]^V 
= V+KiU^^^^V*ulfv, 

(3 3) 

K2 = V* + V*uiyU^2^K2 = V* + V*ull^KiU^l^V* 
= V*+K2UlfvU^2^V*. 
Note that formally one has the power series expansion 

oo 

Ki = Y,{vugW*ullYv, 

oo 

m=0 

which is rigorous if ||[/22'll < 1 or ll^if II < 1- These power series expansions combined with the 
superposition principle were used in [^] to give a formal proof that the composition rule for 
scattering matrices was given by the generalized star product. 

With these preparations the generalized star product U = U'-^^ *v U^^^ of the unitary matrices 
[/'^' and ?7*^^ is defined as follows. Write the («i +n2- 2p) x {ni +n2- 2p) matrix ?7 in a 2 x 2 
block form as 

rr^fUn Un 

\U2l U22 



g 
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where Un is an {rii- p)x{ni- p) matrix, U22 is an («2 ^ p) x ('^2 - p) matrix etc. These matrices 
are now defined as 

Un = ul\> + U[l%uiyui]', 
(3.4) '^^^ 

Ul2-U^^2 f^2Ui2 < 

U2i = U^fK,U^\\ 

In particular for an arbitrary nxn unitary matrix U and all p such that I < p <n the 2px2p 
matrices E = ^ ^ serve as units for the byproduct when V = 1, 

Further we will need the following Perron-Frobenius-type result which for the sake of gener- 
ality will be formulated to cover also the infinite-dimensional case: 

Lemma 3.2. Let a compact operator A on a separable Hilbert space be a contraction, i.e. 
< 1. Suppose that X=\ is an eigenvalue of A. Then 

( i) every c & such that Ac = c also satisfies A*c = c and hence also A* Ac = AA*c = c, 

( ii) the geometric and algebraic multiplicities of the eigenvalue X = I are equal. 



Proof. The claim (i) is an easy consequence of the singular values decomposition (see e.g. [g5 



p. 155]). Thus we have Ker(A - 1) = Ker(A* - 1) = (Ran(A - 1))"^. The claim (ii) now follows 
from the fact that the geometric and algebraic multiplicities of an eigenvalue X are unequal iff 
Ran(A - X) n Ker(A - X) is non-trivial. □ 

Also we will make use of the following 

Lemma 3.3. Let A and B be linear compact operators on a separable Hilbert space ^ such that 
< 1, ||B|| < 1 andABb = bfor some b & f). Then B*Bb = b. 



Proof. Without loss of generality we may assume that \\b\\sj = 1. By Lemma ^ 
(3.5) B*A*ABb = b. 



Therefore by well-known inequalities for the singular values of compact operators (see e.g. [33]) 
we have 

1 < siAB) < siA)s{B) < \\A\\ \\B\\ < 1, 

where s{K) denotes the maximal singular value of a compact operator K, i.e. the maximal eigen- 
value of the self-adjoint non-negative operator K*K. This gives s{AB) = s{A) = s{B) = 1. From 
s{A) = s{B) = 1 it follows that X = 1 is a maximal eigenvalue of A* A and B*B. By the min-max 
principle (see e.g. [ 51 , Theorem XIII.2]) any c g Sj, ||c||f, = 1 maximizing (c, K*Kc) < 1 satisfies 
K*Kc=c. Moreover, if (c, ^r*A'c) = 1 with some ||c||i^ < 1, then ^*A'c=c and ||c||^ = l. Therefore 
since 

{Bb,A*ABb) = {b, B*A*ABb) = 1 

we obtain 

(3.6) A*ABb = Bb 



and ||BZ7|| = 1. The relations (3.5) and (M) imply that B*Bb = B*A*ABb = b. □ 
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Suppose now that the unitary matrix U^^^ is not V-compatible with U^^\ In this case the Unear 
subspaces of C 

e = Ker(I - VU^2^V*ull^), C = Ker(I - U^\^V*ul^^V), 

B = Ker(I - V*u[fvU^\^), % = Ker(I - u[fvU^2^V*) 
are nontrivial. By Lemma 32 we also have 

e = Ker(I - Ulfvu^fv*), C = Ker(I - V*ulfvU^f), 

'B = Ker(I - U^\^'"v*ulf'"v), 'B = Ker(I - VU^fv*ulf*). 

Obviously Q=V*C and 'B=V'B. Since V is unitary this implies dim C=dim C and dim 'B=dim "B. 
Furthermore we have 

Lemma 3.4. The subspaces 23 and C have equal dimensions, dim S = dim C. Moreover 



S = lin span{V*[/*f c, c e C} = lin span{[/^2' ^ g C}, 

' e = lin spa.n{V U^\^b, b & B} = lin span{[/{f Vz^, b&'B}, 

and 

B = hn span{f/*iV?, ?g C} = lin spanjyLf^j^V, J'g C}, 

e = lin span{[/^2'^*^, Z^g S} = lin span{V*[/{f zJg S}, 

f/Zj [/^f c = 0/or a// c g C, 
f/i/j U[2^b = 0/or a// g S, 
(iv) U^fc = Oforallce 6, 
fvj U[l^*b = 0forallb& S. 

Proof. Let c,- g C, / = 1, . . . , ^ < > 1) be a (not necessarily orthogonal) basis in C. For all 
/= 1, . . . ,^ we have 

(3.7) {I-VU^2^V*u[^^)ci = 0. 
Multiplying these equations by V*u[^^ from the left we obtain 

{i-v*ulfvu^l^)v*u[^^Ci = o. 

Thus 

(3.8) lin span{V*[/|f c, c g 6} c S. 
By Lemma it follows from (3^) that 



(3.9) Ulfu[fci = ci, i=l,...,k. 
Hence 

(3.10) dimhnspan{V*C/{fc, cGe} = dime 



and therefore by (3.8) dim C < dim 23 



Let bi G C^, i=l, . . . ,k' < phe some basis in "B. We have 
(3.11) {I-V*ulfvU^l% = 

for all / = 1, . . . ,k'. Multiplying these equations by VU22 we obtain 

(I - VU^2^V*u{f)VU^2^bi = 0, 
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and thus 

(3.12) lmspan{VU^2^b, g S} c 6. 
Again by Lemma 33 it follows from ( |3.1 1| ) that 

(3.13) U^l^*U^2^bi = bi, i=l,... ,k'. 
Thus 

(3.14) dimlin span{VU^2^b, g S} = dim S 



and therefore by (3.12) dim B < dim C. So we have proved that dim 23 = dim C. The inclusion 
( p^ ) and the equality ( p.lOD imply that 

lin span{y*[/{f c, c g 6} = 

The inclusion ( 3.12 ) and the equality ( [3.14| ) imply that 

lin span{V U^l^b, Z> g S} = C. 

The proof of the relations 

S = lmspan{U^\^*V*c, cgC}, 

e = hn span{[/{f VZ7, b e "B} 

is similar and will therefore be omitted. 

We turn to the proof of (ii) - (v). By the unitarity of U^^'> from it follows that J/jf V^f c;= 
0. Since Ker A*A = Ker A for any linear operator A we obtain the claim (ii). By the unitarity of 

from ( pj| ) it follows that U^^ U^n^i = which proves the claim (iii). 
As already noted the vectors c,- and bi also satisfy 

(I - ulf*VU^fv*)ci = 0, (I - U^fv*ulf*V)bi = 0. 



A final application of Lemma |33| yields 

vug%g^*v*ci = a, v*ulfulfvbi = bi 

which by the unitarity of U^^\ U^^^ and V implies (iv) and (v) which completes the proof of the 
lemma. □ 

Lemma 3.5. (i) The matrices I - VU^22^*uf-^ and I - U^^^*VU^^*V* map bijectively onto 
itself, 

(ii) the matrices I - 'V*u'^^'VU22 <^nd I - U22 V map S-*- bijectively onto itself, 

the matrices I - W2^V*U\\'V and I - V*U\\' VU^l^ map C bijectively onto itself, 
(iv) the matrices I — ?7|jV?722^V'* and 1 — ^22 V*U^^ map S"*- bijectively onto itself. 

Proof. Since V is unitary by the definitions of 23 and C it suffices to prove (i) and (ii). By the 
definition of C we have that 

(c, VU^^Vuf^a} = {Ulfvug^*V*c, cj_) = {c, CA_) = 

for any c G C and any c_l g C"^. Thus (I - VU^\^V*uI]\a_ g 6"^ for all c_l g e^. Conversely, 
by Lemma 3^ for any c± g C"*" there is a unique d & Q-^ which satisfies the equation (I - 
VU22V*u[Y)d = c±. This proves the claim (i). The claim (ii) is proved similarly. □ 

Theorem 3.6. If at least one of the off-diagonal blocks of the matrices U^^^ and U^^^ (i.e. U^2' 
U21, U[^, or U21} of maximal rank, then the matrix U^^^ is V -compatible with U^^^ for all 
unitary px p matrices V. 
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Proof. We recall that fov p<n\l2 the {ni-p)xp matrix is not of maximal rank (= min{?ii - 
p,p}) iff there is a vector b & such that ul^b = 0. For p > ni/2 the matrix is not of 
maximal rank iff there is a vector c g such that c = 0. 

Let us suppose that the matrix U'-^^ is not V-compatible with ?7*^'. Then by Lemma 3A it 
follows that all off-diagonal blocks of ?7*^' and U^^^ are not of maximal rank. □ 

Actually we have also the following result. Let a unitary nxn matrix U be written in the 
block form 



U = 



JJ21 U22^ 

where Uu is an (n — p) x{n- p) matrix, U22 is a. px p matrix etc. with p being an arbitrary 
integer such that l< p<n. 

Lemma 3.7. The matrices Uu and U21 are simultaneously either of maximal rank or not of 
maximal rank. 

Proof. Let us suppose that p < nil. Then the {n- p)x p matrix U\2 is not of maximal rank iff 
there is a non-zero vector b & C such that U^b = 0. From the unitarity of U it follows that 

C/i*i Un + 1/21 U22 = 0, [722^/22 + UnUu = I, 

and therefore 

C/2*lf/22^ = 0, U;2U22b = b. 

Thus U22b and U2ib ^ Ker j . From this it follows that the (n- p)x p matrix i/ji is not of 
maximal rank, and thus the px{n- p) matrix U21 is also not of maximal rank. 

Now let us suppose that n> p> nil. Then the matrix U\2 is not of maximal rank iff there is a 
nontrivial vector b g C"~'' such that ?7j*2^ = 0. Again from the unitarity we have 

f/21 t/n + U22UI2 = 0, Un Un + UnU^ = I, 

and therefore 

U2iUnb = 0, UnUnb = b. 
Thus, the px{n- p) matrix U21 is not of maximal rank. □ 

We will show now that the ^-product can be extended to arbitrary, not necessarily V-compatible 
unitary matrices. We will prove that the operators 

and 

u['2'K2 = - v*u[yu^^r'v* 



are actually well-defined. From Lemma ^ it follows (see [ )36| , Section 1.5.3]) that 

Pqi{I - VU^lW*ulfr^ and P^i{I - V'VlfvU^l^T^ 



are well-defined. From (ii) and (iii) of Lemma 'iA it follows that C c Ker U^^ and c Ker UI2 
and thus 

u^]\i - vu^l^v*ulf)-^ = u^]^Pqi{i - vu^2^v*ulfr^ 

and 

ulWi - v*uiyu^\Y = - v'v[yu^\Y 

are well-defined. Similarly one can show that the operators 

U^fvil - C/^2'^*C/{i V)-! and uH^il - ulfvU^l^y^ 
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are also well-defined. With this we obtain that the relations (34) indeed also define the general 



ized star product of two non-compatible unitary matrices. Moreover, we have 
(3.15) 



f/^i v(i - [7^2 v^ii c = 0, ui^'Vii - ui^^'vu^^ v r b = ^ 

for all c G e, G c G e, and ^ g S. 

Theorem 3.8. For arbitrary unitary matrices U''^\ U^^\ and V the matrix U = *y ?7^^' is 
unitary. 



This theorem was proved in Appendix C of [42] in the case when is V-compatible with 



U^^K For the general case the proof is given in Appendix A below. 

Analogously one can prove associativity of the generalized star product. More precisely let 
?7*^' be a unitary ni, x «3 and V' a unitary p' x p' matrix with p' < n2, p' <ni,. If p + p' <ni, then 

t/'^' *v iU^'^ *v' U^'^) = *v U^'^) *v' U^'^ 

holds. 

Theorem 3.9. The generalized star product is a continuous operation in each of its two argu- 
ments, i.e. for any unitary matrices U'-^\ U^^\ U^^\ and V such that U^^^ and U'-^^ have equal 
size there is a constant C > depending on U^^^ and V only such that 

where the norm \\ ■ \\ is an arbitrary matrix norm. A similar estimate holds with respect to the 
first argument. 



In [ ]42| ] we proved that the scattering matrix of a self-adjoint Laplacian on an arbitrary graph 
is a continuous function of X > 0. Theorem together with the composition rule given in 
Section Q below allows to give an alternative proof of this fact. We will not dwell on the details 
here. 



In the sequel we will use the notion of the Moore-Penrose pseudoinverse (see e.g. [|70|]). 
Recall that for any (not necessarily square) matrix M its pseudoinverse M* is uniquely defined 
by the Penrose equations 

MM*M = M, M*MM* = M'', 
{M*Mf = M*M, (MM*)* = MM*. 

One also has 

M** = M*\ 
RanM* = RanM*, 
KerM'' = KerM*, 

and MM* = PRanM, M*M = PrsuM', where denotes the orthogonal projection onto the linear 
subspace Ji. Moreover 0* = 0. If M is a square matrix of maximal rank then M* = M~^. 

Let U again be an arbitrary unitary nxn matrix written in the block form with some l<p<n. 

Lemma 3.10. IfKerUn = {0} then Ker {U21 - U22U^JUuy = {0}. 

Proof. Assume the converse, i.e. let there be c g C^, c^O such that 

iU2i-U22U^2UnTc = 0, 

or, equivalently, 

U;,c-U*nU^2*U22C = 0. 
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We multiply this equation by U21 from the left and use the unitarity of U which in particular 
implies 

t/2lt/21 + t/22t/22 = I- 

This yields 

(3.16) C - U22U22C - U2lU^iU^2*U22C = 0. 

Again by unitarity we have U2iUyi = -?/22^^i*2- Recall that 

UnUu = {UnUnT = d - PK.rUnf = I 
by the hypothesis of the lemma. Thus, from ( |3.16[ ) it follows that c = 0. □ 

Now we turn to a discussion of the inverse of a unitary 2px2p matrix U with respect to the 
generalized star product *p, i.e. the existence of the unitary matrices and such that 

U^*p U = U*p U'^ = E, 

where E is the 2px2p matrix ^ ^ (in the pxp block notation). We will not discuss general 

necessary and sufficient conditions for the existence of and U^, but simply restrict ourselves 
to a special case. We will prove 

Theorem 3.11. Let U be an arbitrary 2p x2p (p> \) unitary matrix. Let at least one of the 
pXp matrices U\2 and U21 be of maximal rank {=p). Then there exists a unique unitary 2p X 2p 
matrix U' such that 



(3.17) 



U' *pU = U *p U' = E. 



(3.18) 



Proof. By Lemma both matrices U12 and U21 have maximal rank. We will discuss only the 
second of the relations ( 3.17 ). In block notation this relation has the form 

Un + Uu{I-U'nU22r'u'uU2i=0, 

U^2 + ^2l(I - U22U'n)-'U22U[2 = 0, 
Uu{I-U'nU22r'U[2 = I, 

U^,{I-U22U'nr'U2i = I. 

By Theorem 3.6 the matrix U must be compatible with U' such that (I - U[^U22)~^ and (I - 
U22U[i)~^ are both well defined. We multiply the first of the relations ( 3.18 ) by U^2 from the 
left. Next we multiply the resulting expression by I - U[^U22 thus obtaining 

U'n{U2i - U22U^l Un) = -U^l U, i . 

By Lemma 3.10 we have that U2\ - U22Lfx2^n is invertible and thus 

(3.19) U[, =-U^jUu{U2i - U22U^jUur'. 



From the third relation in (3.18) we obtain 

U[2 = {I-U[iU22)Ui^ = Uij + [/ri[/ii([/2i - U22UijUnr'U22Uil 



The fourth relation in ( |3.18| ) gives 

(3.20) t/^i = U2l{I - U22U'n) = (f/21 - U22U{^Unr'. 



The second relation in (3.18) determines 1/22- 

It remains to prove that U' is unitary. By the unitarity of the matrix U we have 

(C/2*i - [/nt/iVV2*2)(t/2i - U22U^2Un) = 1 + t/nt/iV't/il ^11- 
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Therefore 



([/2*i - UuUn V2*2r^[I + t/nt/r2 'f/i1t/ii](t/2i - U22UilUnr^=l 



and thus by ( p9| ) and ( |3^ ) we obtain 

u[,*u[,+u^;u^,=i. 

The relations U[2*U[2 + ^^22*^22 = ^^11*^^12 + ^^21*^22 = 0' and U[2*Un + ^^22*^21 = can be 
proved similarly. 

The left inverse is constructed similarly and by means of the obvious relation 

{U21 - U22U^lUnr'U22U^l = t/2"/t/22(t/21 " U22U^^Uur'U22U^j 

is easily shown to be equal to the right inverse. □ 

Corollary 3.12. Let G be the set of all 2p x 2p unitary matrices with px p blocks U\2 and U21 
both being of maximal rank, endowed with the generalized star product *p as a multiplication 
and E as a unit. Then G is a group isomorphic to U(/j, p). 

The proof will follow from the arguments given in Section ^. In particular, the group iso- 
morphism between G and U(/7, p) is given by the formulas and ( 5.13| ) below. We note that 



this isomorphism generalizes the well-known set isomorphism between the group SU(1, 1) and 
a subgroup of SU(2). 

The set of all 2p x 2p unitary matrices endowed with the generalized star product *p as a 
multiplication and E as a unit is no longer a group but only a semigroup. 

4. Composition Rule for the Scattering Matrices 

Now we apply the generalized star product to prove the composition rule for the scattering 
matrices on graphs. For this we only need the special case V = Ip, the px p unit matrix, and 
so we introduce the notation *p:=*v=ip- Let Fi and F2 be two graphs with ni> I and ?i2 ^ 1 
external lines, respectively, labeled by £1 and E2, i.e. #(£1) = ni, #(£2) = «2 and an arbitrary 
number of internal lines with given fixed lengths (see Fig. |l|). Furthermore at all vertices we 
have local boundary conditions giving Laplace operators A(ri) on Fi and A(F2) on F2 and the 
scattering matrices SiCk) and S2{^)- Let now £j and £2 be subsets of £1 and £2 respectively 
having an equal number l<p< min{?ii , ^2} of elements. Also let cpo: £? — > £2 be a one-to-one 
map. Finally to each k & £j we associate a number a^ > 0. With these data we can now form a 
graph F by connecting the external line e £ j with the line (^aik) g £2 to form a line of length 
aic. In other words any interval [Ok,°°k), ^ g £? belonging to Fi and the interval [0^,g(k),°°(?o{k)) 
belonging to F2 is replaced by the finite interval [Ok,ak] with 0/t being associated to the same 
vertex in Fi as previously and a^ being associated to the same vertex in F2 as Otp^^k) before in 
the sense of the discussion at the end of the previous section. Recall that the graphs need not be 
planar. Thus F has n = n\+n2-2p external lines indexed by elements in (£1 \ £'j') U (£2 \ £2) 
and p internal lines indexed by elements in £° in addition to those of Fi and of F2. We denote 
this set by J12 such that the set of all internal lines of the graph F is given by J = Ji U J2 U J12. 
There are no new vertices in addition to those of Fi and F2 so the boundary conditions on Fi 
and F2 define boundary conditions on F resulting in a Laplace operator A(F). Suppose that the 
indices of £j in £1 come after the indices in £1 \ £j (in an arbitrary but fixed order) (see (3.1 )). 
Via the map cpo we may identify £2 with £'j' so let these indices now come first in £2, but again 
in the same order. Finally, let the diagonal n2 x n2 matrix V{a) be given as 

, / exp 0\ 
^^^=[ " ij' 

where exp(/\/Xa) again is the diagonal px p matrix given by the p new lengths au, k g £3. 



SCATTERING MATRICES ON GRAPHS 



15 



Recall that by Theorem 2A all eigenfunctions of the operator -A(r) have the form 

fo, 7g£, 
\aje'^' + ^je-'^\ j&3, 

where the coefficients aj and satisfy the homogeneous equation 

ZA.B,ai^) (^aj =0, a = {a,}™ 1 , p = {^jY^U 

with the matrix Za^b.oO^) defined by (^). We define the linear subspace HuO^) of C""'"^'" as a 
set of all vectors (0, a, P)^ for which aj = = for all j g 3 12, 

(4.1) = {( = {0, a, p)^ G Ker Za,b,JX) c C"+''" | ay = p,- V; e On} ■ 

Obviously for X t. OA,B.a this subspace is trivial, i.e. iini'k) = {0}. Let T(r, J12) c R be the 
set of those eigenvalues of -A(r) for which Ker Za^b.oO^) © ''^12(^) is nontrivial. Obviously the 
eigenfunctions corresponding to the eigenvalues from T(r, J12) have nontrivial overlap with J12, 
i.e. supp(|) n J12 has non-zero measure. 

Let S(ri,r2) c R+ be the set of those X > for which KeriV{g)S^22('^)V{q)Sf^{X) - 1) is 
nontrivial. 

Theorem 4.1. With the above notations E(ri, r2) = T(r, 3 12) D M+. The composition rule 

(4.2) S(k) = SdX) *p V{a)S2{XWia) 

holds for all X g M+. IfX g T(r, J12) and X>0 then its multiplicity equals 

dim Ker(-A(ri )-X) + dim Ker(-A(r2) - X) 
+ dimKer(V(a)42 (>^)V'(£)5fj'(>i) - 1). 
In particular if the eigenvalue X>0 is such that X g T(r, J12) then its multiplicity equals 

dim Ker(-A(ri )-X) + dim Ker(-A(r2) - X). 

Here dimKer(-A(r/) - X) denotes the multiplicity of the eigenvalue X regardless whether it 
is embedded into the absolutely continuous spectrum or not. 

Note that by Lemma ^ dim Ker{V{q)S22{X)V{q)Sf^{X)-l) equals the algebraic multiphcity 
of the eigenvalue ^= 1 of V(£)4y(?t)V(£)5fi'(?i). 

If by cutting p internal lines of an arbitrary graph F with local boundary conditions, the 
graph will be decomposed into two disjoint subgraphs Fi and F2, by Theorem ^[T] the scattering 
matrix Sr can be obtained from the scattering matrices and ^r, at the same energy. Thus, 



using (4.2) iteratively the scattering matrix associated to any graph can be obtained from the 
scattering matrices associated to its subgraphs each having one vertex only. In fact, pick one 
vertex and choose all the internal lines connecting to all other vertices. This leads to two graphs 



and the rule (4.2) may be applied. Iterating this procedure L times, where L is the number of 



vertices, gives the desired result. 



Proof of Theorem 4.1. We split the proof into several steps. 

1. First we suppose that Si{X) is compatible with V {q)S2{X)V (a) and prove that the composi- 
tion rule (U) holds. Let \|/J(x, X; F/), j g £/ U J/ for any k& £/ be the solution of the Schrodinger 



equation with the operator -A(F/), / = 1, 2 (see (|2.4|)) at energy X. Let ^i{x, X) be «/ x«/ matrices 

(4.3) ['¥i{x,X)]ji^ = \ii){x,X;ri), j,k&Ei, 1 = 1,2 
such that 

(4.4) »P/(x, X) = e-'^^'I + e'^^'SiiX), 1=1,2. 
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Fig. 1 . Gluing of two graphs. 



Observe that and are linearly independent functions and therefore the scattering 

matrices may uniquely be recovered from ^i(x, 7i). The columns of ^i(x, 7i) define the external 
part of solutions of the Schrodinger equation for -A(r/) at energy X.. We are looking for a square 
matrix 

such that its {n\ - p) + («2 - p) = ni + n2 - 2p columns defines the external parts of a solution 
to the Schrodinger equation for -A(r). Here the indices are assumed to be arranged such that 
the first indices are those of £i \ followed by the indices of £2 \ £2- The aim is to obtain 
^(x, X.) from ^\(x, A,), ^2(^, X), and the lengths a = {as}sG3i2 of the new internal lines J12. By 
the above observation this will determine the scattering matrix S(k). The strategy will be to find 
new solutions of the Schrodinger equations for -A(r/) with incoming plane waves m 
the channels £/ \ £? which agree suitably in the channels k& 8-^ and (po(fc) e £2. 
With the conventions made above we write 

(^{ni-p)x(ni-p) 

1 5^"^"^^^^(X.) 



(4.5) 




where the superscripts denote the sizes of the blocks. For arbitrary px{ni- p) matrices C\ and 
C2, respectively, consider the n\ x {n\ - p) and n2 x (ni - p) matrices 



(4.6) 



^2{x,X;C2) = e-'^ (^^^ +e'^S2{X) (^^^ 



Here I stands for the {ni- p)x{ni - p) unit matrix and stands for the (?i2 - x ("1 - p) zero 
matrix. Obviously, the columns of ^lix, X; Ci) are the external parts of linear combinations of 
the columns of ^/(x, X; F/), and thus define the external parts of solutions of the Schrodinger 
equations for the operators -A(r;), 1=1,2. Note that 5>i(x, X\ C{) has an incoming plane wave 
in any of the channels k& £1 \ £? and ^2{x, X; C2) has no incoming plane wave in all channels 
kG £2\£^. 
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Now we make the coordinate transformations x ^ ak-x on the Unes (po(^)£ £2 ^i)- "^^^ 
reason for this is as follows. Under the gluing process Fi , r2 — > F (see Fig. [l]) the two half-lines 
corresponding to ^ g £j and cpo(^) g £2 ^•'^ replaced by the interval [0, a^-], giving the new lines 
in J12. This may be realized by identifying a point P on the half-line corresponding to k & £*' 
and with coordinate x{0<x<ak) with the point Q on the half-line corresponding to cpo(fc) g £2 
with coordinate - x. In particular x = akon the ^-line corresponds to x = on the cpo(^)-line 
and vice versa. Applying this transformation to <I>2(x, X; C2) we obtain in this new coordinate 
system 



^f{x,X;C2) = 



i\/X(a^x) 






i\/%X.a—x) 



1 



We now require that <I>i(x, X;Ci) and <I>2-'(x, X; C2) agree on the lines labeled by J12. This will 
fix C\ and C2. Indeed, we then obtain 



iv/^.cPX("i-/')(X) + ^.'v^^f (X)Ci 



^C2 + e 



iVX(a-x) qpxp 



By the linear independence of the functions e'^" and e it follows that 
(4.7) 



Ci=e'^^i5f'' (X)C2, 



and thus 



(4.8) 



C2 = 



I-f''^^5f''(X)e'^^^5f^(X) 



Ci=e'^^5f'' (X) 



I - e'^^f (?i)e'^^5f ^(X) 



(X). 



Since for any invertible A and U the identities ?7A ^ = (A?7 ^) ^ and A = {U ^A) ^ hold, we 
have 



e'^^5f''(X)e'^^i5f^''(X) 



= e 



e 



i\f%a qP^P 



I-5f^(X)e'^5f''(?i)e^ 



^/Xa 



Since S\(k) is assumed to be compatible with V {q)S2(k)V {q) the inverses in (18) are well de- 
fined. 

Similar to ( p3| ) and according to the ordering convention made above we write the scattering 
matrix S{X) for the graph F in the block form 



S(k) = 



^(n 1 -p)x{ji 1 -p) Q^-^ ^(n 1 -p)x{;i2 -p) Q^^ 
5(n2-p)x{"i-p)(^) 5'{'i2-p)x{;i2-p)(^) I ' 



where the superscripts denote again the sizes of the blocks. Since <I>i(x, X; C\) has an incoming 



plane wave in any of the first «i channels k& £1 \ £*', equations (4.8) allow one to determine 
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^(«i-/;)x(«i-p)(^^ and 5("i-/')x(«2-p)(X): 



_ ^(ni-p)x(n,-p) 



5' 



,py.(ni-p) 

1 



(n2-p)xp Q^^^i^/}ia 



= S 



jpx(;ii-p) 



To determine the remaining blocks of the scattering matrix S{X) instead of ( [4.6[ ) we consider 
the ni X (?i2 - and ni x (?i2 - p) matrices 



^i{x,X;Ci) = e 



(4.9) 



Ci 



<I>2(;c, X; C2) = e-'^' (^^ ] + e'^'^52(X) 





Ci 

C2 



with arbitrary px{n2- p) matrices C\ and C2. Again 0/(x, X; C/) are the external parts of some 
solutions of the Schrodinger equations with the operators -A(r/), 1= 1,2. Now ^\{x, X; Ci) has 
no incoming plane waves in the channels ^ g £1 \ g^, but <I>2(.x, X;C2) has an incoming plane 
wave in any of the channels g £2 \ £2- 

Repeating the arguments used above we obtain the following matching conditions for Ci and 
C2: 

(4.10) 

and thus 



C2=e'^^5r'ci, 



5f ^(X)C2 + SP^'-P'iX) = e-'^^^Ci, 



pxini-p)/^ X _ „-i-/kqj 



Ci = 



I - e'^5f ^(;i)e'^5f ^(X) 



jV)ia(;PX{n2-p) 



C2 = e'^^5f"(X) [l-e'^-^5f''(X)£''^^5f^(X) 



a), 

i\fXac;py-(n2-p) 



Since Si(k) is compatible with V {d)S2QC)V {q) the inverses are again well defined. From this and 
from (0^) it follows that 



^(n2-p)An2-p)Q^^ = 

= s\ 



in2-p)An2-p)Q^^ + 4"2-P)^''(X)C2 
in2-p)An2-p)Q^^ + 4"2-'"^^(X)e'^^5f^''(?l) 



I - f^'^-^^f ''(?l)e'^^5f ^(X) 



i\/XaQP'x(n2-p) 



^(ni-p)x(«2-/?)^^^ _ 



_ Aiit-pjxp 



I - e'^^f (X)e'^5f ^(X) 



i\/Xaf;PX(n2-p) 



By the definition of the generalized star product (3.4) we obtain (4.2). 
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2. Now suppose that X g S(ri, r2). We prove that the composition rule remains valid. 
Also X G T(r, J12) and the multipicity of X equals 

dim Ker(-A(ri + dim Ker(-A(r2) - X) 

(4.11) +dimKer(y(a)42'W(a)5fi'(?i)- 1). 

The assumption X g S(ri, r2) implies that 



and 



have nontrivial kernels. This implies that the homogeneous form of the equations ( |4.7| ) and 

Ci=e'^5f^(X)C2, 



(4.12) 
and 
(4.13) 



5f^(X)Ci=£'-'^C2, 



C2 = e'^5f^(X)Ci, 



5f''(X)C2 = e-'^^Ci, 



respectively, have nontrivial solutions. It is easy to prove that the inhomogeneous equations 
( p^ ) and ( 4.10| ) still have solutions in this case. Consider for instance the equation (4/7), which 
is equivalent to 



C2 = e 



{X) + e'^^S':''''{X)e'^S'^''{X)C2. 



By the Fredholm alternative this equation has a non-trivial solution iff 

(4.14) 5f<'"""\X)*e-'^Z7 = 
for any O^b & C satisfying 

5f ^(X)*e-'^Sf ^(?i)*e-'^^Z^ = b. 
By Lemma ^ with A = 5f ''''(X)* and B = g-'^^S'^'''' (X)* e''"^-!- we have 

(4. 15) e''^SP''PiX)SfP{X)*e-'^b = b. 
From the unitarity of V{g)Si{X)V{g), which states in particular that 



e'^Sf''(X)5f'' (X)*e-'^^^ 



I, 



and from ( 4.15| ) it follows that 



J \/%a ^px{n I -p) Q^^P^'^n 1 -p) 



\Xte 



b = 0. 



Since Ker C*C = KerC for any operator C we obtain ( 4.14 ). Equation ( 4.1C ) is discussed simi- 
larly. 



From ( |4.6| ) and ( |4.9| ) it follows that the Schrodinger equation with the operator -A(r) for 
given value of the spectral parameter X>0 has (nonunique) solutions which have the form 



(£>iix,X;Ci) = e 



(4.16) 



Ci 



+ e 



'\/X>c 



Sx{X) 



^2{x,X;C2) = e-'^ (^'- 



+ e 



SiCX) 





Ci 

C2 
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and 



^i{x,X;Ci) = e-' 



(4.17) 



+ e'^^\Si(X) 



(£>2{x,X;C2) = e 



-i-/Lx ( C2 




+ e 



'\/Ajc 



S2{X) 



Ci 

C2 





where Ci and C2 (Ci and C2) solve ( [4.1 2D (( [4.13[ ), respectively). Note that Ci = Ci and C2 = C2. 
On the lines in the set J12 the quantity <I>i coincides with <I>2 and and 02- We will now prove 
that 



(4.18) 



,{ni-p)xp 



{X)Ci=0, 



An2-p)xp 



(k)C2 = 0. 



Thus, the functions ( |4.16 ) and ( 4.17 ) are zero on all external lines of the graph T and their 
support has nontrivial overlap with the interval lines J12. 

From ( PI ) it follows that Ci = e'^5f''''(>i)e'^^^5f (X)Ci. By Lemma U we have 



By unitarity it follows that 



5f^(X)*5f"(X)Ci=Ci 



and thus S 



,{ni-p)xp 



(k)C[ = 0. The second relation in (4.18) is proved similarly. 



Now we note that from (4.12) and (4.13) it follows that 



RankCi = RankC2 = dim Ker ( £''^^i5f''(?i)f''^^^5f' '(X) - 1 



The columns of ( [4. 16[ ) correspond to linear independent eigenfunctions of -A(r) for the eigen 



value X. There are precisely dim Ker {^e'^^S2^'^{X)e'^^S^^^{X) - ij such eigenfunctions and 
the supports of all them have nontrivial overlap with the internal lines J 12. 
3. Let X G T(r, Jn) n M+ and let 

(4.19) dim {KerZA,B.a(k) Q UniX)) = k, 

where the linear subspace LuiX) is defined by (p~T|). This means that there are precisely k 
eigenfunctions of -A(r) which disappear if we cut the internal lines Jn- We will prove that 
Xg sen, r2) and that 

(4.20) dim Ker(V(a)4y(?t)V(a)5fi'(>i) -l) = k 

which in turn implies that 

dim Ker(-A(r) -X) = dim Ker(-A(ri )-X) + dim Ker(-A(r2) - X) 

+ dimKer(y(£)42(>^)V'(£)5fi'(?i) - 1). 

From the existence of the above mentioned eigenfunctions it follows that these eigenfunc- 
tions can be constructed by means of superposition and matching of the solutions (4^) of the 
Schrodinger equation for the operators -A(ri) and -A(r2) at energy X> 0. For any vectors 
Ci, C2 G the functions 



^i{x, X,Ci) = e 



+ e 



SiiX) 



^2{^> X,C2) = e 



_ -iVAj: / C2 







+ e 



S2{X) 
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define the external parts of solutions of the Schrodinger equations for the operators -A(r/), 



/ = 1,2. Since the eigenfunctions are supported on internal lines of the graph F (Theorem |2.1[ ) 
the vectors Ci and C2 must satisfy 



j{ni-p)xp 



such that (|)i(x, X, Ci) vanishes in any of the channels k & £1 \ £j and ^2ix,)^,C2) vanishes 
in all channels g £2 \ £2- Making the coordinate transformation x —^ a - x on the lines 
(po(^) e £2 s £j) and requiring that (|)i(x, X, Ci) and (|)2(£- x, X, C2) agree on the lines labeled 
by J 12, we obtain 



Ci=£''^5f^(?l)C2, 

C2, 



5f ^(?i)C- - --'^^ 



^2 
: e 



or equivalently 



(4.21) 



e'^5f "(X)e'^^5f ^(X)Ci =Ci, 
C2 = e'^5f^(X)Ci. 



Linear independent solutions of ( 121 ) correspond to linear independent eigenfunctions of -A(r) 
and vise versa. Thus the condition ( 4. 19| ) implies (120). This completes the proof of the theo- 
rem. □ 



Note that if F is simply the disjoint union of Fi and F2, i.e. if no connections are made 
(corresponding to p = and n = ni + n2), then S{X) is just the direct sum of Si{X) and 5'2(A,). 
Also V*S(k)V = S2^^{X) *v S(k) for any scattering matrix with n open ends and any unitary nXn 
matrix V, where 



I 

1 



Similarly SiX) *v S^i^Ck.) = VS{X)V* 



Fig. 2. The graph from Example | 



Example 4.2. Consider an arbitrary self-adjoint Laplacian A(A, B) with local boundary con- 
ditions on the graph depicted in Fig. H where the distance between the two vertices is a. The 



composition rule (4.2) with 



V{a) = 



22 



V. KOSTRYKIN AND R. SCHRADER 



easily gives 



(4.22) 



Sn 


_ c(l) 


+ s 


S22 


_ c(2) 
"^^22 


+ s 


Su 


_ c(l) 




S21 


_ c(2) 
-"^21 


4?( 



M) M) 2iaVls.-l 

^i22^ne ) , 

(I) M] 2ia^/I^-l 



" ^^22^^! 



)" 



where the S-matrices are written in the form analogous to (3.1) 



- S\ 



(1) 



^21 



/c.(2) r.(2)' 
c(2) _ "^ll ^12 
- I j,(2) ^(2) 
^^22 , 



-"21 



leaving out the X— dependence. These relations are equivalent to the well-known factorization 
formula [57, 58, 47|, 64, 55] applied to the Laplacian on a line with boundary conditions 

posed at x = and x = a. 




Fig. 3. The graph from Example p3| . The arrows show the positive direction 
for every edge. The edges 3 and 6 have the length a and the edges 4 and 5 the 
length b. 



Example 4.3. Consider the graph depicted in Fig. ^ where the length of the edges 3 and 6 
equals a and the length of the edges 4 and 5 equals b. Let the boundary conditions be given as 

¥i(0) = ¥3(0) = V4(0), 

¥2(0) = \|/3(a) = ¥5(0), 

^^'^^^ ¥'i(0) + ¥'3(0) + ¥'4(0) = 0, 

\|/2(0) + V5(0)-V3(«) = o. 

\^4W = ¥6(0) = 1^7(0), 

¥5(^) = ¥6(«) = ¥8(0), 
-V4(Z7)+V6(0) + V7(0) = 0, 
-V5(ft)-V6(a)+X|/8(0) = 0. 

Obviously they define a self-adjoint operator which we denote by A(a, b). The scattering matrix 
corresponding to this operator (as defined by (13) and ([2.4Dj will be denoted by Sa,h(X)- To 
determine this 4x4 matrix we first consider the graph depicted in Fig. ^ where the length of 
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the edge 3 is supposed to be equal a. The boundary conditions ( 4.23| ) determine the self-adjoint 
operator. The corresponding scattering matrix we denote by Sa(k). From (2.7) it follows that 



Fig. 4. The graph from Example [^.3| . The arrows show the positive direction 
for every edge. 



(4.24) 



— ( J^/T^ 



9e 



-4 e 
2{e'^' - 3e-'^) -4 
V -4 2(e'^ - Sf--'^^ 



+ 3e 



l{e 



By Theorem 4. 1 the scattering matrix Sa^hQ^) Is given by 

(4.25) Sa,b(^) = SaiX) *2 V{b)Sa{'k)V{b), 

where 

V{b) = diag(e'^^ 1, 1), b = {b, b) g M^. 



We now compute the 2x2 matrices K\ and K2 entering the definition (3.4) of the generalized 
star product, thus obtaining 

^-1 = = (^e'^' - 9e-'^'^ L, 

(L)i 1 = (L)22 = e^'^il - e^'^) + 9e-^'^{9 - e^'^^) - 2(9 + Ue^'^), 
(L)i2 = (L)2i = Se^'^^'ie'^' + 3e-'^'). 
From this it follows that 

det Ky^ = det = (^e'^ - Qe"'^ 

.g-4,VI« ^(^^2 _ g4)(^ _ g)2 ^ i6t^(_256 - 128^ + 44^^ _ 3^3-, 

where t, = exp{2iVXa} — 1 and T\ = exp{2iVXb} — 1. Obviously these determinants vanish 
jj:- ^2i\/Xa _ ^2i\/Xb _ j ^^^^ show that there are no other zeros. Note that the embedded 
eigenvalues of the operator —A{a, b) are determined by the equation g2i\/Ia _g2(\/X/5 _ ^ iji^i 
for incommensurable a and b there are no embedded eigenvalues. 
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For e^iVXa = g2iVXfc _ | ^j^^ matrix Sa(k) is not compatible with V{b)Sa(k)V{b) and 

where +1 corresponds to exp{/-\/Xa}=+l. Obviously Ker A'j"^ =Ker A'j^ is the subspace spanned 
by the vector (1, =Fl)^- Further, 

(«.),.(;)4(> -)(>,)=o 

Thus, as proved in Section ^ the generalized star product is well defined also in the case when 
the matrix Sa(k) is not compatible with V (b)SaO^)V {b). 




Fig. 5. The graph from Example p3[ The arrows show the positive direction 
for every edge. The edges 3 and 6 have the length a and the edges 4 and 5 the 
length b. 



Example 4.4. Consider the graph depicted in Fig. ^ where the length of the edges 3 and 6 
equals a and the length of the edges 4 and 5 equals b. Let the boundary conditions be given by 

¥i(0) = ¥2(0) = >|/3(0), 
¥4(0) = \|/5(0) = \j/3(a), 
Vi(0) + V2(0) + V3(0) = 0, 

V4(0) + V5(0)-V3(a) = 0, 
¥4(^) = ¥5(^) = ¥6(0), 
\|/6(a)=\|/7(0) = \|/8(0), 
-V4(^)-V5(ft) + V6(0) = 0, 

-<(a) + V7(0) + i|/8(0) = 0. 

Obviously they define a self-adjoint operator which we denote by A(a, b). The scattering matrix 
corresponding to this operator (as defined by (2.3) and (2.4)j will be denoted by Sa,hO^)- To 



determine this Ay. A matrix we first consider the graph depicted in Fig. ^ where the length of 
the edge 3 is supposed to be equal a. The boundary conditions (4.23) determine the self-adjoint 



operator. The corresponding scattering matrix Sa(^) can be obtained from (124) by means of 
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Fig. 6. The same graph as in Fig. Hbut with different ordering of the external lines. 



the permutation of its lines and columns thus giving 



Sa(k) = ( e'^ - 9e-' 



2{e 



i\f%a 



-4 
-4 



+ 3e 

-4 

-4 



-4 
-4 



2{e 



3e- 



-4 

_ -4 

) e' 



By Theorem 4.1 the scattering matrix Sa^hO^) i^ given by 

(4.26) SaA'^) = SaCk) *2 V{b)Sa{'k)V{b), 

where 

V{b) = diag(e'^^ e'"^'', 1, 1), b = {b, b). 



We now compute the 2x2 matrices K\ and K2 entering the definition (3.4) of the generalized 
star product, thus obtaining 

= = l^e'^' - 9e-'^y^ L, 

(L)ii = (L)22 = e^'^il - 5e^'^^) + 9e-^'^^{9 - 5e^'^) - 18(1 - e^'^), 
(L)i2 = (L)2i = -4e^'^{e^'^' - 9e-^'^). 
From this it follows that 



det K^^ = det K2 



1 _ fj^i 



+9^e-^'^{9 - lOe^'^^ + e'*'^^) - 36e^'^^{\ - 6e^'^'' + 5e^'^') 



_jg2g-2,-VX«(9 _ i4g2,Vw> + 5e^'^) + 18(27 ■ 



Obviously these determinants vanish if e^'^' = 1. One can show that there are no other zeros. 
Note that the embedded eigenvalues of the operator —A{a, b) are determined by the equation 

For e^i^A^ = I f/jg matrices Sa(k) and V{b)Sa(kWit) '^f^ not compatible and 



KV=K:^ =-4 



(giVXfl _ 9g-i\/Xa^2 \^1 \J ' 
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Obviously Ker ^Tj ' = KerKj ^ is the subspace spanned by the vector (1,-1)^. Further, 



and 



{Samn \ .\)=-4(j I) (_\)=0 



{V{b)Sa(k)V{b))n [ _\) =-Ae^'^' (\ J) (\ ] =0. 



Thus, as proved in Section ^ the generalized star product is well defined also in the case when 
the matrices Sa(k) and V{b)Sa(k)V{b) are not compatible. 



As already discussed in p2[ ] multiple application of (4.2) to an arbitrary graph allows one 
by complete induction on the number of vertices to calculate its scattering matrix from the 
scattering matrices corresponding to single- vertex graphs. If these single vertex graphs contain 



no tadpoles, i.e. internal lines starting and ending at the same vertex, then ( |4.2| ) give a complete 
explicit construction of the scattering matrix in terms of the scattering matrices for single vertex 
graphs. In case when a resulting single-vertex graph contains tadpoles we proceed as follows. 
Let the graph T have one vertex, n external lines and m tadpoles of lengths a, . To calculate the 
scattering matrix of T we insert an extra vertex on each of the internal lines (for definiteness, 
say, at X = a,72). At these new vertices we impose trivial boundary conditions corresponding 
to continuous differentiability at this point. With these new vertices we may now repeat our 
previous procedure. Thus in the end we arrive at graphs with one vertex only and no tadpoles. 

5. Special Case n\=n2 = 2p: Transfer Matrices 

This section is devoted to the construction of the transfer matrix for Schrodinger operators 
on graphs with an even number of external lines. The transfer matrix formalism for general 



Schrodinger operators on the line is well known (see e.g. [19]). Its relation to the scattering 
matrix is discussed in e.g. [ |7T| , pT] ]. 



Fig. 7. The graph with n = 2 and m = 0. 

We start with the simplest example of a Laplace operator on the graph with n = 2 and m = 
(see Fig. ^ which is equivalent to a Schrodinger operator on the hne with point interaction. The 
boundary conditions given by the relation 



(5.1) 

where the matrix 



(^2{0)\^j,(a b\ (^\fdO)\ 
[^^'^{0)J \c d) ' 



a b 
c d 



SL(2,M), 

and fj is real, lead to self-adjoint Laplacians (see [^, 68|, |16|, \[% |3||). Conversely, from the 



viewpoint of the von Neumann extension theory (see e.g. [|59|]) relation ( |5.1[ ) describes almost 
all (with respect to the Haar measure on U(2)) self-adjoint Laplacians A{A,B). If exp{2/^} = 1 
the operator A(A, B) is real, i.e. commutes with complex conjugation. In particular, the choice 
a - I = d - I = b = 0, exp{2//7} = 1 corresponds to the 5-potential of strength c (see e.g. [Q]). 
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By definition tlie transfer matrix is a 2 x 2 matrix M(k) g U(1) X SL(2, R) satisfying 
(5.2) 



;l.n^f¥l(0)^_^¥2(0)^ 



In fact, in the case at hand it is given explicitly as follows 

If exp{2//j} = 1 the matrix M{X) is unimodular, i.e. M{X) g SL(2; M). 

The transfer matrix possesses the following equivalent description. Any solution of the 
Schrodinger equation with the operator -A(A, B) for the energy X > has the form 



From this and (jS^) it follows that there is a matrix A{X) g U(1) x SU(1, 1) c U(l) x SL(2; ( 
(with the inclusion in the group-theoretical sense) such that 



(5.3) 
and 



ay 



-iVl 



1 



-iVl 



For X>0 the matrix A(X) is related to the scattering matrix 



^^'''^^>-\T2{X) L{X)) 



by the relation 



I 1 m \ 

T2{1) T2(k) 

L(k) mX) 

\T2iX) \T2(kt) 



where 



T,(k) 
T2(X) 
R{X) 
L{X) 



2e'^{a - ibVX + ic/VX + d)'^ , 

le-'^'ia -ibVX + ic/VX + dy^ , 

{a - ibVX + ic/VX + d)~^ {a -ibVX- iciVX - d), 

{a -ibVX + ic/VX + d)~^ {-a - ibVX - ic/VX + d). 

Note that T[{X) = T2{X) for all X > if the operator A{A,B) is real, i.e. exp{2//7} = 1. This is in 
analogy with Schrodinger operators on the line with potentials which are necessarily real (see 
e.g. giU]). 

The factorization rule from Example ^ can now be written in the form 
(5.4) A{X) = A^^\X)U{a)A^^\X)U{a)-^ , 

where 

/ „-i\f\a 



U{a) = 








The relation (5^) is the special case of the well-known factorization formula ^ |58|, ^ 
64 , applied to the Laplacian on a line with point interaction. 
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It is easy to realize that the transfer matrix cannot be defined for arbitrary boundary condi- 
tions. For instance, the Dirichlet (\|/2(0+) = \|/i(0+) = 0) or Neuman (\|/2(0+) = \\i\{0+) = 0) or 
mixed (\|/2(0+)+fc2¥2(0+)=Vi(0+)+^i\|/j(0+)=0) boundary conditions introduce the decouphng 
A(A, B) = Ai ©A2, where Aj, j = 1,2 are the Laplacians on L^(0, 00) with corresponding boundary 
conditions. Recall, however, that the scattering matrix is well defined even in these cases. The 



composition rule ( [4.22D (see Example ft.2D remains valid. 

Now we consider an arbitrary graph F with an even number of external lines n = 2p. We 
enumerate the external lines in an arbitrary but fixed order. The external part of an arbitrary 
solution of the Schrodinger equation with -A(A, B) at the energy X > has the form 



(5.5) 

We define the transfer matrix 



(5.6) 



uj{x) = aje''^ + i 



,n. 



\bj 



/Vi\ 



\ an j 



To prove that A(X) is correctly defined it suffices to show that for arbitrary constants {aj,bj), 
7 = 1,... ,p there is a solution to the Schrodinger equation with the operator -A(A, B) whose ex- 



ternal part has the form ( p.5[ ) and this solution is unique up to its internal part. The external part 
of any solution to the Schrodinger equation with the operator -A(A, B) is a linear combination 
of the columns of the matrix-valued function 



(5.7) 



A{X) 



Thus, the columns of ( |5.7| ) have to satisfy (p.q), i.e. 

'SuiD 5i2(X)\ / I 

J \S2l{X) 522(?l) 

where the px p block notation is adopted. Writing A{X) as 

■^W=fA"Al a' ml 

V^A2l(X) A22(X)y 

we obtain 

Au(X)5ii(X) + Ai2(X) = 0, 
Aii(X)5i2(X) = I, 

A2l(^)5ll(X) + A22(X) = 52l(X), 
A2l(^)5l2(^) = 522(X). 

Let us suppose that det SuC^) ^ 0. Then 

Au(X) = 5i2(X)-\ Ai2(X) = -5i2(X)-'5ii(X), 

A2l(^) = S22{l)SM-\ A22(X) = 521 (X) - 522(X)5i2(X)-l5ii(X). 

Thus, we proved that for det 5i2(?t) ^ the transfer matrix exists and has the form 



(5.8) 



(5.9) 



>12W Oil 
52l(^)-522(X)5i2(X)-l5ii(X) 



Also, its definition (| 
(5.10) 



immediately leads to the following factorization formula 
A(;i) = A^^\X)U{a)A^^\X)U{ar^ , 
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where the diagonal unitary matrix U (a) is given by 



U{q) = 







Note that formal arguments based on the superposition principle leading to ( 5.10 ) have appeared 
earlier in [20|. As for related results we mention that in [50] it was shown that the transfer matrix 
of a Schrodinger operator on the line with a matrix-valued potential can be written in the form 



3- 

Lemma 5.1. //det Sn(k) * then K(k) g \J{p, p). 



Proof. Obviously the coefficients a 1, 


. . . ,an,bi,.. 












\ =S(k) 






ytn) 





From the unitaiity of the scattering matrix it follows that 

|aip + . . . + |a„p = \bi\'^ + . . . + [ 

or, equivalently, 

\ai\^ + ... + \ap\^ - \bi\^ - \bp\^ = \bp+i\^ + . . . + 
This relation and ( ^ ) complete the proof of the lemma. 
Let us summarize the above results of the present Section: 



\ap+i I 



\a„\ 



□ 



Theorem 5.2. If det 5'i2(^) ^ then the transfer matrix A(k) g \J{p,p) as given by (5.9) exists 
such that for an arbitrary {ai, . . . ,ap,b\, . . . ,bp) g C" there is a solution of the Schrodinger 
equation with — A(A,B) at the energy X > whose external part has the form (5.5) and the co- 
efficients {bp+i, . . . ,bn,ap+\, . . . ,an) g C" are given by (5.6). The composition rule for the 
scattering matrices (4.2) is equivalent to the multiplication formula ( 5.10| )/Qr the transfer ma- 
trices. 

In addition for real operators we have 

Theorem 5.3. If the operator A(A,B,a) is real and if in addition det S\20^) ^ then A(?i) g 
SU(A7;>)cSL(2p;C). 



Proof. From the well-known determinant formula for block matrices (see e.g. [30, Section II.5]) 



(5.11) 



det 



'An Ai2 
^^21 A22 

which follows from the decomposition 



= det All det (A22-A2iAi}A 



12) 



All 
A21 



A12 
A22 



it follows that 



By Theorem I 



All 
A21 



det A(k) = 



I 



A7}A 



11^12 



A22-A2lAi}Ai2 

det 521 (X) 



det Si2iX)' 

I we have SuCk)^ = SziCk) and thus det A{X) = 1. 



□ 



We turn now to a discussion of the assumption det Si2(k)^0. For the scattering matrix of the 
graph depicted in Fig. ^ with the boundary conditions ( |4.23| ) (see Example O) det (Sa(X))j2=0 
for all X>0. 
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Theorem 5.4, Suppose that det Si2(k) = 0. Then the transfer matrix A.QC) exists such that for 
arbitrary 



(ai, . . . . . . ,Z7p) G Ran ( !| ^ ) 



there is a solution of the Schrddinger equation with —A{A,B,a) at energy X>0 whose external 
part has the form ( |5.5| ) and the coefficients {bp+i, . . . , bn, ap+[, . . . , a„) g C" are given by (5.6). 

Proof. The external part of any solution to the Schrodinger equation with the operator -A(A, B, q) 
satisfying the conditions of the theorem is a linear combination of the columns of the matrix- 
valued function 

1 



(5.12) ^{x,X) 



P, 



{Ker5'i2{?^)) 



where *F(x, X) is given by (5^). Thus, the columns of ( 5.12 ) have to satisfy (p^), i.e. 

'AniX) Au{X)\ (SniX) SniX)P^K.rSna)A = ( ^ ^(Ker5„a))^ 
^A2iiX) A22{X)) \ I y V^2l(^) S22(X)P^^,,Sn(^)K 

The solution of this equation can be wiritten in the form 

Aii(^) = Sn{X)\ An(X) = -Sn(XfSn(X), 
A2i(?t) = S22(X)Sn(X)\ A22 = 521 (X) - S22{X)Su{XrSn{X), 
where ★ stands for the Penrose-Moore pseudoinverse. □ 

Note that any vector of the form (c, 0)^ with c g KerSi2(X) satisfies A(X) = 0. Thus 
det A{X) = 0. 

Inspection of the proof of Theorem 5.4 shows that the transfer matrix cannot be extended to 
a subspace larger than 

1 

^(KerSi2{?.))^, 

If det A'^^\X) = det A^^\X) = then Ran U {a)A^^\X)U (a) and Ker A(i)(X) may have a nontrivial 
overlap and therefore the multiplication formula ( |5.10| ) does not hold in this case. 

Example 5.5. Consider the graph depicted in Fig. ^ with the boundary conditions from Exam- 
ple |4. j[ For all X G M+ such that e^"^ = \ we have that Ker S\2{X) is nontrivial and 

1/1 



Ran 



■fKer5i2(>t) 2 I zpl 1 

p 



where +1 corresponds to exp{/\/Xa} = ±1. Suppose that 



{a\,a2,b\,b2)^ g Ran 



P{}isrSn(k))^ 



{a\,a2, b[,b2)^ g Ran 



or, equivalently, 

'0 

yO PKerSnW 
In particular, we can choose 

ai=a2 = 0, b[ = \, b2 = ^\. 

It is easy to check that there is no solution to the Schrddinger equation with these boundary 
conditions. 
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Example 5.6. Consider the graph depicted in Fig. g with the boundary conditions as in Exam- 
ple 4.4 KerSi2(X) is nontrivial for all X g M+ and 

1 / 1 -A _ 1 



-1 1 



(KerSi2(X))-L 



Suppose again that 



and choose 



{ai,a2,bi,b2) i Ran 



a\ = a2 = 0, 



= 1, 





°(KerSi2(>.))-l 



b2 = \. 



Again it is easy to check that there is no solution to the Schrddinger equation with these bound- 
ary conditions. 



The statement converse to Theorem 5.2 immediately follows from Theorem 5.4 



Theorem 5.7. If the transfer matrix A(k) exists in the sense of Theorem 5^ then det SuO^) ^ 
and the corresponding scattering matrix is given by 

-Anil)-' Anil) AuHr' 



(5.13) 



S{X) = 



A22(^) - A2i(^)Aii(^)-'Ai2(^) A2i(X)A„(X)- 



Proof. Suppose that det S^O^) = 0. Then by Theorem we get det A(k) = 0, which is a 
contradiction. Thus, det Si2{'^)^0 and therefore by Theorem 52 det A{'k)j^O. The representation 
Kl^ follows from (K§). □ 



Appendix A 

Here we give the proof of Theorem 3.6 which claims that for arbitrary unitary matrices U^^\ 
U^^\ and V the matrix U = U^^^ *v U^^^ defined by (3^) is unitary. As already noted in [42] it 
suffices to prove only the relations 

t/i 1*^/11 + ^/21*^/21 = I, 

(A.l) 



The remaining relations 
(A.2) 



Uu U12 + U21 U22 = 0. 
Ui2*Un + U22*U22 = I, 

Uu*Un + U22*U2i=0. 



follow immediately from ( |A.1[ ). To see this for an arbitrary unitary matrix U we define an 
involutive map U ^ given as 



U = 



Un U12 



U22 U21 



JJ21 U22J \Ul2 Uii 

Direct calculations show that the following "transposition law" 

(A.3) [/- = [/(2)%^,[/a)^ 



holds whenever U = U^^^ *v U^^\ Assume that (AJ) holds for arbitrary unitary U. Replacing 
the matrix U by given by ( |A.3| ) transforms the relations ( |A. 1[ ) into (A.2). 

By the definition of the generalized star product (3.4) and by the unitarity of U^^^ the first of 



the relations ( |A. l] ) is equivalent to 
(A.4) 



U21 U21 +U21 V ^^2^12 ^11 +^11 t^l2 ^2tVii Vt/ji 
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Since the opposite case was already considered in [ |4^ ] we further assume that the matrix U^^^ is 



not V-compatible with ?7*^'. From Theorem 3^ it follows that all off-diagonal blocks u[2, ^2^, 

u[^, and 1/2^ are not of maximal rank and thus Ker U21 is nontrivial. Let di, l</<^=dim Ker U21 

be an arbitrary basis in Ker U21. From the unitarity of the matrix ?7*'' we get UI2 U^^dj = for 
all l<i<k. Thus 

"^21 ^21 +^21 ^ ^11 ^2t^l2 ^11 +^11 ^12^2^/11^1/21 



di = 



for all 1 < / < /c. Hence to prove ( |A.4| ) it remains to show that 
(A.5) 



^21 ^21 +^21 ^ ^11 ^2^2 ^11 +^11 t/12 ^2t/ii V/t/21 



d = 



for any d g ( Ker U21 ) = Ran 1/21 ■ Therefore we set d = U21 d, where J g is an arbitrary 



vector. Thus, the relation (|A.5[ ) holds whenever 
(A.6) 



^21 ^21^21 +^21 ^^11 ^2^2 ^11^21 +^11 ^712^21/11^^/21^/21 



+ ^21 ^^11 ''^2^12 ^12 ^2t/ii V/t/21 ^21 +^21 ^1^21 U21K1U21U21 



d = 



for all d G CP. 



First we note that by Lemma 3.4 (iv) the relation (A.6) holds for all J g C. Therefore it 



suffices to prove that ( |A.6D holds for all d g S"*-. Observe that in this case by Lemma 3^ (i) and 
by the unitarity of the matrices ?7*^' and ?7*^' we have 

u[fu[\Xfd = -UlfulMTd = -U^*d + t/l^lf Jg S\ 
(A.7) J= V J- Jg %\ 

U^2iU2i *d = d-U^^U^2Td& e^. 



To prove the first relation in ( A/7) it suffices to show that for any J g C"*" and any b & H 
(A.8) -{b, U^22*d) + {b, U^l^*U^\^U^fd) = -{U^2^b, d) + {u!22U^2*^22^' d) = 0. 



I — I n^* A 1 

By the definition of S and by Lemma 6^22 ^^22^ ~ ^ ^"^y t> ^ H which proves ( |A.8| ). To 



prove the second relation in ( [4.7[ ) it suffices to show for any J g C"*" and any b & "B 
(A.9) 



{b, u[fVd) - {b, ul^^VU^^U^fd) = 



iv*ull^b, d) - {u!2-^u^^ V*i7{f b, d) = 0. 



By Lemma (i) V*uf^b g C. By the definition of C and by Lemma |3| [7^2^ [7^2^* ?= c for any 



c G C which proves ( [4.9[ ). This also proves that 



(C, J) - (C, [7^2'^^r^) = (C' ^) - (^22^^22**^, d) = 



for all J G C and all c g C from which the third relation in {k.l) follows. 
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From Lemma ^ and the definition (3^) of the matrices Ki and K2 it follows that 

(i) Ki maps C"*" onto C"*" bijectively, 

(ii) ^j* maps C"*" onto C"*" bijectively, 

(iii) K2 maps S"*" onto S"*" bijectively, 

(iv) K2 maps 23"*" onto 'B'^ bijectively. 



(A. 10) 



Noting that ul^/u^^b^ = b±- U^l^*U^2^b^ g "B^ and U^f*U^fci_ = c_l - [/{f V{f ci g 6-^ 
due to (kJ_) and ( A.IO ) we can write the l.h.s. of (A^) in the form 

U^f[-I- V*u[fK2U^f - U^l^K2U{y + V*ulfK;K2U[y 



Similar to [42| one can easily prove that for any b± g S-*- and c± g C"*- the following relations 
hold: 

K*iKi^_L = Cj_ + [/^2'^2f/i*i + V*ulf*K2U^f^_i_ 
+V*u{fK*2 U^fu^^K2U[y^^, 

Klulfull^KiCi, = V*a + ulfK;U^22*V*)Ulfulf{I + VU^l^K2Ulf)V^^, 

K*2K2bl_ =^j_ + /:2^22^V*[/{f * + [/{iV[/^2'^2^-L 
,*,,(1)*,.*,,(2)%,(2) 



V*U\\' Ul\'VU^^'K2b^. 



Inserting these relations in ( |A. 1 1| ) with the choice c± = U^^d and b_\_ = U\\'Vci_ we obtain 
that it vanishes thus completing the proof of the first relation in (A.l). The proof of the second 



relation in ( A. 1 ) is similar and will therefore be omitted. 
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